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Abstract 

Experiments on tiie cuprate superconductors demonstrate that these materi- 
als may be viewed as a stack of Josephson junctions along the direction normal 
to the Cu02 planes (the c-axis). In this paper, we present a model which de- 
scribes this intrinsic Josephson coupling in terms of incoherent quasiparticle 
hopping along the c-axis arising from wave- function overlap, impurity-assisted 
hopping, and boson-assisted hopping. We use this model to compute the mag- 
nitude and temperature T dependence of the resulting Josephson critical cur- 
rent jc{T) for s- and d-wave superconductors. Contrary to other approaches, 
d-wave pairing in this model is compatible with an intrinsic Josephson effect 
at all hole concentrations and leads to jdT) oc T at low T. By parameteriz- 
ing our theory with c-axis resistivity data from YBa2Cu307_5 (YBCO), we 
estimate jc{T) for optimally doped and underdoped members of this family. 
jc{T) can be measured either directly or indirectly through microwave pen- 
etration depth experiments, and current measurements on Bi2Sr2CaCu208 
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and La2-xSrxCu04 are found to be consistent with s-wave pairing and the 
dominance of assisted hopping processes. The situation in YBCO is still un- 
clear, but our estimates suggest that further experiments on this compound 
would be of great help in elucidating the validity of our model in general and 

the pairing symmetry in particular. 

PACS numbers: 74.50.+r, 74.20.Mn, 74.25.Nf, 74.25.Fy 
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I. INTRODUCTION 



While much of the early work on the cuprates has focused on the copper oxide planes, 
recent experimentalS and theoreticallHi efforts are beginning to address the c-axis properties. 
In large part, this interest is motivated both by the availability of large, high-quality single 
crystals and by the realization that a full understanding of the superconductivity cannot be 
reached until the role of the third direction is established. Indeed, the third dimension is 
ssential for the high critical temperatures in some theories. i0 

While a complete theoretical description of the c-axis properties of the cuprates is clearly 
desirable, a full understanding of these properties is hampered by their unusual experimental 
signatures. In the normal state, for example, the resistivity along the c-axis (pc) is quite 
generally metallic in the optimally doped materials and becomes semiconducting when the 
hole doping is reduced, while the in-plane resistivity (pab) always shows metallic behavior.lll 
At first sight, the metallic behavior of both pab and pc in the optimally doped compounds is 
understandable within the conventional Bloch-Boltzmann theory of transport; however, the 
semiconducting temperature dependence of pc with metallic pab in the underdoped samples 
is difficult to reconcile with this picture. In the superconducting state, also, measurements of 
hysteretic current-voltage curves, Fraunhofer-like c-axis critical currents as a function of in- 
plane magnetic field, and emission of microwave radiation at ac Josephson effect frequencies 
strongly indicate that the cuprates can be viewed as a stack of SIS Josephson junctions 
along the c-axis.ini This situation is not at all consistent with the properties of conventional 
superconductors. 

One picture which may account for these anomalous properties is based on incoherent 
c-axis transport. In conventional metals, the presence of the crystal lattice leads to the for- 
mation of three-dimensional Bloch waves. The basic postulate of the incoherent approach is 
that Bloch waves in the c-direction do not form due to the layered structure, intrinsic dis- 
order, and/or strong electronic correlations in the cuprates. This situation is not equivalent 
to localization of the quasiparticles in the c-direction, though, because they may still hop 
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incoherently between the Cu02 layers. Working under this assumption, we will develop a 
theory which explains the unusual temperature-dependence of pc in the normal state and 
the intrinsic Josephson effect in the superconducting state and which may also shed light on 
the question of the order parameter symmetry in these systems. 

The specific model we adopt was first articulated by Rojo and Levini and asserts that the 
Cu02 layers can be thought of as independent 2D Fermi liquids weakly coupled along the c- 
axis by three processes: direct inter-layer quasiparticle hopping,!'! hopping assisted through 
static disorder ,0 and hopping assisted through dynamic, boson-mediated (e.g., phonon- 
mediated) scattering.!!! The semiconducting temperature dependence of pc is then associated 
with the reduction of the dynamic inter-planar scattering as the temperature is lowered. The 
c-axis conduction is therefore analogous to that of a very dirty, but nevertheless metallic, 
system.! 

We will show that all three of these processes give rise to Josephson coupling, although 
their relative importance depends on the material and on the order parameter symmetry. 
It should be no surprise that all three mechanisms produce Josephson coupling between the 
layers; direct and impurity-assisted hopping processes formally resemble the tunneling pro- 
cesses considered in conventional SIS junctions,&^i and boson-assisted hopping has been 
known to contribute to tunneling in superconducting junctions in both the quasiparticle0 
and Josephson0 channels for some time. In the cuprates, interest in these mechanisms has 
been piqued by the interpretation of experiments involving the c-axis transport! and flux 
quantization in superconducting rings. This interest is reflected in a rapidly developing 
theoretical literature. In particular, we mention related results on intrinsic Josephson cou- 
pling in the cuprates by Graf et a/., El which came to our attention after completion of this 
work. As in this paper, these authors employ an incoherent hopping model to describe the 
inter-planar coupling in the cuprates. We differ in that Graf et al. consider only impurity- 
assisted hopping while we also include direct and boson-assisted processes. As a result, we 
find a non-zero (i-wave Josephson critical current in fully oxygenated YBCO without the 
need to introduce an anisotropic impurity scattering rate. In addition, because the present 



approach can describe the temperature dependence of the normal-state transport, we have 
been able to compute the full temperature dependence of the critical currents. 

We present our calculations as follows. In the next Section, we describe our model 
and the approximations used in our calculations. For later reference, we will also discuss 
the normal-state c-axis conductivity cXc and derive expressions for the contributions to this 
conductivity by the three processes. We then move on to the superconducting state and 
derive the general forms of the Josephson critical current jc which arise from each inter-layer 
hopping mechanism. The direct and impurity-assisted critical currents are found to have 
familiar forms, but the boson-assisted contribution gives rise to a novel expression. Finally, 
we quantify our model by parameterizing experimental resistivity data in the YBa2Cu307_5 
(YBCO) system and using the estimated parameters to predict jc{T) for different oxygen 
stoichiometrics and pairing symmetries. Since jc can be measured directly0 or inferred 
from microwave penetration depth experiments,^ our predictions may be compared to ex- 
periment. Based on these calculations, we conclude that low-temperature measurements of 
jc{T) should be able to distinguish s- from (i-wave pairing regardless of oxygen stoichiome- 
try. Moreover, we find that the c-axis penetration depth ratio should have either a 3D BCS 
character if the pairing is rf-wave or a 3D BCS character which crosses over as a function 
of hole doping to a 2D Ambegaokar-Baratoff character if the pairing is s-wave. Although 
experiments have not converged in the YBCO system, we discuss the implications of our 
model for the pairing symmetry in the other cuprates in the last subsection. A preliminary 
report of these results has been published elsewhere.ill 



II. INCOHERENT HOPPING MODEL 

A. Description and Assumptions 

We begin by writing the Hamiltonian for the electronic system in terms of the Hamil- 
tonians of the individual Cu02 layers Hm and an inter-layer coupling Hamiltonian H±: 
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Hc\ = J2m Hm + Hj^. Here, the term "layer" denotes a composite object consisting of all the 
Cu02 planes within a single unit cell, so the transfer Hamiltonian H± takes quasiparticles 
from one unit cell to another. Motivated by Ref. ^, we take H± = A + A\ with 

ima j 

and Cjmo- a quasiparticle annihilation operator for site i in layer m and spin projection 
cr. The terms in H± represent a direct hopping arising from wave-function overlap {t±), 
an impurity-assisted hopping due to static disorder of strength Vim, and a boson-assisted 
hopping represented by the field (pim with coupling strength Qim- The total Hamiltonian is 
then the electronic Hamiltonian added to the Hamiltonian governing the dynamics of the 
field: i/tot = i^ei + H^. 

In studying the properties of this model, we make two assumptions; the first relates to 
how we include incoherence in the model and the second to how we view the intra-layer 
dynamics. Microscopically, incoherence may arise from strong intra-layer scattering,i a 
non- Fermi-liquid ground state within the layers,! dynamic inter-layer scattering,i or strong 
electromagnetic fluctuations.! We do not attempt to construct a fully microscopic theory 
of this incoherence, but simulate it by taking iJ^ to be a small perturbation on Hm and 
calculating its effect within second order perturbation theory. This approach is similar in 
spirit to that of loffe et a/.il and should yield the correct features of a more complete theory. 
In order to perform quantitative calculations within this model, we also assume that the 
normal state of Hm for each layer m is a Fermi liquid and that the superconducting state 
may be described within Migdal-Eliashberg theory.!! 

From these two assumptions, the procedures outlined in the rest of this paper may be 
applied to any type of inter- layer disorder Vim or inelastic scattering (pim- For concreteness, 
however, we specify these parameters. First, we take the disorder to be delta-function 



correlated: Vim = and VimVjn = V^dijdmn, where the overbar denotes the usual average over 
impurities.@ Moreover, we take the electron-inter-layer-boson coupling to be structureless 
(fi'jm = 9) and the bosons themselves to be represented by an Einstein spectrum of frequency 
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Qq. These choices do not affect the character of our results in s-wave superconductors, but 
they do affect the critical currents in ci-wave superconductors in a way which will be discussed 
in Section III. 

When it is necessary to include intra-layer scattering processes in our calculations, we do 
so through a temperature-dependent scattering rate h/Tab(T) which is linear in T at optimal 
doping, at least down to the critical temperature Tc. In addition, all wave-vector integrals 
are performed by the usual technique of restricting the wave vectors to the Fermi surface, 
taking the density of states to be constant, and then integrating over the remaining energy 
dependence. This approximation is appropriate if the Fermi level is far from any structure 
in the density of states (e.g.. Van Hove singularities) and the other energy scales in the 
problem are much less than the band width. We assume these conditions hold here. 



B. Normal-State Conductivity 

To illustrate the application of our model and for future reference, we compute in this 
subsection the normal-state c-axis conductivity within the incoherent hopping model out- 
lined above. By assumption, this calculation is to be carried out to second order in H± and 
thus becomes formally identical to computing the conductivity of an NIN tunnel junction, 
with H± playing the role of the transfer Hamiltonian.iii Applying the standard theory of 
tunneling,lli we compute the correlation function X(r) = — (Xr A{t)A''{0) ^ with A defined 
by Eq. (|l]) and with the r-dependence of A determined by Y.m Hm + H^f,. The diagrams used 
in this calculation are shown in Fig. |l|; note that vertex corrections to X{t) are of higher 
order in the inter-layer hopping amplitudes and so are neglected. Analytically continuing the 
Fourier transform of X{t) to real frequencies gives the c-axis current density produced by 
a voltage V across nearest neighbor layers as jciV) = — (2e//ia^) lmX{eV), from which the 
conductivity cTc = c {djc{V) / dV)\v=o is obtained. (Throughout this paper, a (c) denotes the 
intra- (inter-) layer unit cell dimension). The resulting conductivity is the sum of three terms, 
each corresponding to one of the inter-layer hopping processes: cTc = a'^^^'^^* + cr™P + cr™*'^. 
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The first term in cJc results from the direct hopping characterized by t± and can be 
written 

ar-^ = aoN,tl g), (2) 

where Nq is the density of states per unit cell per spin at the Fermi surface and ao = 
{ATce"^ /h){c/a'^), Up to factors of order unity, this conductivity reproduces the results of 
Kumar and Jayannavari and loffe et a/.ii Note that a'^^^'^^* is proportional to the intra-layer 
lifetime Tab, just as in the usual Drude expression for the conductivity, despite the fact 
that the conductivity is viewed as a tunneling process. To understand this surprising result, 
recall that the Kubo formula for the conductivity in 3D metals leads to a particle-hole bubble 
with group velocity factors at the vertices, while these vertices become the tunneling matrix 
element Tk'k for tunneling calculations. The Kubo conductivity thus has a different form 
than the tunneling conductivity in 3D materials. For the 2D metals comprising the layers in 
our model, however, the tunneling matrix element for direct hopping is T^'k = tj_6k'k, and so 
vertices are diagonal in 2D wave vector. Consequently, the Kubo and tunneling formalisms 
give the same results.^ 

The second and third terms in ac result from the assisted processes and are found to be 

a'r = cro{NoVy (3) 

for impurity-assisted hopping and 

for boson-assisted hopping. Observe that the impurity-assisted hopping conductivity is the 
usual NIN tunneling conductivity.El Also note that the procedure in Ref. ^ of simulating 
the presence of boson-assisted hopping in (Xc by allowing the impurity-assisted hopping 
amplitude to become a function of temperature is microscopically justified by the more 
complete calculations which yield Eqs. 

The form of the conductivity given by Eqs. (^-(^ provides a simple explanation for the 
anomalous normal-state resistivity in the cuprates. At optimal stoichiometry, we expect the 
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layers to be strongly coupled and nearly coherent, implying that the direct hopping should 
dominate cXc. Consequently, pc oc pab as in conventional Bloch-Boltzmann transport.0 Off 
optimal stoichiometry, the direct inter-layer coupling may weaken relative to the assisted 
processes, making the material increasingly incoherent. While this change would not strongly 
affect Pab, Pc would now be dominated by the assisted hopping, the inelastic component of 
which would freeze out at low T [cf. Eq. (^)]. Thus, as T decreases, the net inter-layer 
hopping amplitude would be reduced and pc would increase, as observed experimentally. 
These qualitative statements can be made quantitative through detailed fits of the normal- 
state resistivity to the Rojo-Levin model. 

III. JOSEPHSON CRITICAL CURRENT 

The incoherent hopping model described above views the inter-layer transport as an 
NIN tunneling process. In the superconducting state, this formulation of the transport 
immediately produces an intrinsic Josephson effect between the layers. In this Section, we 
compute the temperature-dependent critical current associated with each inter-layer hopping 
process to determine the general features of the total jc{T). In the next Section, we will 
apply these theoretical results to the cuprates and see what can be revealed about the order 
parameter symmetry and the relative importance of the inter-layer hopping mechanisms. 

A. Formalism 

The technique for calculating the Josephson critical current in an SIS junction is well 
known and is closely related to the corresponding NIN calculation. One simply uses 
the same diagrams from the NIN computation and replaces the normal Green's func- 
tions with anomalous or Gor'kov Green's functions.lii^ii Specifically, we define the corre- 
lation function $(r) = — {T^ [A(r)y4(0)]) and evaluate it using the diagrams in Fig. |1|. 
The Fourier transform of the resulting expression is then analytically continued to real 
frequencies to obtain the Josephson current density produced by an applied voltage V: 
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jj{t) = — (2e/^a^)Im e"^**^^*/^ $(eV^) . Setting = and adjusting the relative phase of 
the order parameters on different layers to maximize jj(0), we obtain the Josephson critical 
current as the sum of three components: 

Jc(T) = jf^^^\T) + fr{T) + jr\T). (5) 

In principle, these critical currents depend not only on the inter-layer hopping mechanism 
but also on the order parameter symmetry.^ 

Although the Josephson critical current and the conductivity are computed from the same 
set of diagrams, there are significant differences. The most obvious is that the propagators 
used in computing the diagrams are different; specifically, the phases of the order parameters 
enter into the Josephson calculation in a fundamental way and affect how the imaginary part 
of the correlation function $ is taken. Moreover, the conductivity is the derivative of the 
correlation function X while the Josephson current is the correlation function $ itself. We 
therefore cannot expect similar behavior from these two quantities. In particular, we will 
see that the boson-assisted hopping does not contribute to the conductivity at low T but 
does contribute to the Josephson current. 

In evaluating the correlation function $, we should ideally solve the full Eliashberg 
equationJni for the layered system and use the results of that calculation to construct the 
anomalous propagators present in $. This approach would naturally include all the effects 
of elastic and inelastic scattering which would naively have a large impact on the Josephson 
coupling. As we will argue below, however, this procedure is not necessary, since inter-layer 
scattering effects are of higher order in the inter-layer hopping amplitudes and the intra-layer 
scattering may be incorporated into a BCS-like gap function. As we shall see, this approach 
is a natural consequence of our model and is supported by the long history of experimental 
and theoretical work on the Josephson effect in SIS junctions. 00 

In performing the reduction from the full Eliashberg equations, we first note that all 
quantities in our model are by construction computed to second order in the inter-layer 
hopping amplitudes. This approximation is reasonable in the superconducting state as long 
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as the hopping amphtudes are smaller than the maximum of the gap function, as we expect 
them to be. Since the vertices in the diagrams for $ are each first order in the hopping 
amplitudes [cf. Fig. 0], the complete second-order calculation involves propagators for only 
single, isolated layers. In particular, there are no vertex corrections to $ and any loss 
of Josephson coupling due to inelastic inter-layer hopping comes in at higher order in the 
inter-layer hopping amplitudes. 

Solving for the propagator is facilitated by observing that the critical current depends 
only on the weak- coupling propagator with the strong-coupling gap function A(k, iujn) (= 
(pCk, iun) / Z (k, iun) in the Nambu notationil) in the approximation where wave vectors are 
restricted to the Fermi surface and the energy integrals are extended to infinity. Hence, all 
intra-layer elastic and inelastic scattering can be incorporated into our model by computing 
their effect on the gap function. Actually performing this computation for the cuprates is 
complicated by the presence of the unusual intra-layer inelastic scattering, which presum- 
ably gives rise to the linear-in-T resistivity. Since the origin of this scattering is presently 
controversial, we are forced to take a phenomeno logical approach at this point. 

Empirically, the cuprates exhibit the generic features of BCS superconductors, although 
with quantitative differences in, for example, 2A/ ksTc. We therefore adopt the conventional 
picture and compute jc using the BCS form for the gap function. Anderson's theorem,ii 
in combination with the experimental fact that even very dirty tunnel junctions exhibit 
Ambegaokar-Baratoff behavior, suggests that this approach should be reasonable in the s- 
wave case. For rf-wave superconductors, the effect of disorder and scattering is more subtle,0 
so our results should only be taken qualitatively. 

We conclude this subsection by recapitulating some of the assumptions which underlie the 
calculations which follow. First and foremost, we have simulated the effects of incoherent 
inter-layer transport through second-order perturbation theory. This treatment does not 
address the details of how the incoherence arises in the first place, so further work will 
be required in order to obtain a fully self-consistent description. We have also neglected 
the detailed structure of the cuprates within a unit cell (multiple Cu02 planes, CuO chains, 
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etc.) and have treated the intra-layer scattering mechanisms approximately; in particular, we 
have neglected the possibility of anisotropic scattering, which may be important in rf-wave 
superconductors. Additionally, we have not attempted a fully self-consistent description 
of the superconducting state, but have relied on BCS theory. As we have argued, these 
assumptions are well justified. We therefore expect that the results which follow will not be 
strongly modified in more sophisticated treatments. 

B. Results 

From Eq. (1), the total Josephson critical current is the sum of three terms, each one 
corresponding to an inter-layer hopping process. Following the preceding discussion, we cal- 
culate in this subsection the critical current from each inter-layer process using the diagrams 
in Fig. and a BCS gap function.!! 

For direct hopping, this procedure gives 

where = + A^, f(x) is the Fermi distribution function, and the angle brackets 
denote an average over the Fermi surfaced In the limit of zero-temperature, j^^^'^^^{0) = 
{2e/ha?) NQt]_ for either s- or (i-wave symmetry, independently oi the magnitude of Ao(0). By 
contrast, jc(0) oc Ao(0) in the Ambegaokar-Baratoff formula.^ As a function of temperature, 
jdirect^ji^^jdircct^g^ varics cxactly as the in-plane BCS superfluid density, which is plotted in 
Fig. 1^ for both s- and rf-wave pairing. At first sight, it is surprising that an incoherent 
treatment of the direct interlayer hopping yields the same temperature dependence as the 
(coherent) in-plane transport. However, as in the normal state, the reason behind this 
behavior derives from the fact that the vertices in the tunneling diagram are diagonal in 2D 
wave vector and so give an effectively 3D result (cf. the discussion after Eq. (|2D). We also 
point out that j^''''^'^*(T) depends strongly on the pairing symmetry at low T, being nearly 
independent of T for an isotropic s-wave gap but proportional to T for a clean (i-wave gap, 
as might be expected from the presence of nodes in Ak. 
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tanh(Ek/2A;BT) ^9/(Ek) 



2Eu 



dEu 



(6) 



Impurity-assisted inter-layer hopping gives rise to a contribution to the critical current 
which has the Ambegaokar-Baratoff form if the order parameter is s-wave,El 

,™.(r).il(.iV„l.)^A„(r)tanh(Mp) (7) 

[cf. the solid line in Fig. 0, and vanishes by symmetry if the order parameter is (i-wave 
and the impurity scattering is isotropic. The assumption of isotropic scattering is crucial 
if j™P(T) is to vanish for (i-wave superconductors; if the scattering is anisotropic, then the 
critical current is in general finite.El One generally expects impurity scattering to be mainly 
isotropic, so j™P(T) is probably small in d-wave superconductors, but detailed calculations of 
the scattering matrix element in these systems are required to make a quantitative estimate. 

Finally, the boson-assisted hopping processes also contribute to the Josephson critical 
current. Direct calculation of the relevant diagram yields 

HP 

jr' = j-,NoXAl{0)lUT), (8) 



where {h = ks = 1) 
'"■^'^ ^ TT Am Jo E Jo 



2 Ag(T) r^o de f'^o de' 
^ A§(0) Jo 'E Jo 'E' 

xl[f{-E')-m] 



+ [f{E') - fiE)] 



{E' + Ef - 91 



{E' + E) coth ( ^"l - coth + ^ 



{E' - EY - 91 



2TJ \ 2T 

fno\ ^ , fE'-E' 



{E' - E) coth {^) -Qq coth 



\2Tj \ 2T 

Qq is the Einstein phonon frequency, A = 2Nog'^ /Qq, and s-wave pairing is assumed. For d- 
wave pairing and structureless electron-boson coupling, the boson-assisted Josephson current 
vanishes .ii 

The T and Qq dependence of /inci(r) are illustrated in Fig. ||for a 100 K superconductor. 
This figure shows several interesting features of li^eiiT). First, we see that this function is 
approximately a non-zero constant at low T, and the magnitude of this constant grows with 
Qq. Second, Iine\(T) resembles the Ambegaokar-Baratoff result for Qq > Ao(0), but exhibits 
a peak when Qq < Ao(0). Finally, near T^., /inei(r) shows a linear dependence on — T. 
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(9) 



We may explain these features simply. First, the Josephson current does not vanish at 
zero temperature because the zero-point energy of the bosons allows the creation of virtual 
bosons to mediate the inter-layer Cooper pair tunneling. Since the zero-point energy is 
proportional to Qq, the increase of /mci(O) with Qq is expected. Second, when Qq > Ao(0), 
Ao(0) is the lowest energy scale and controls the integral, so lir^diT) is structureless. On the 
other hand, when Qq < Ao(0), we see from Eq. (P) that the Qq coth(i7o/2T) terms dominate 
the integral and lead to an increase in Jinci(T) with increasing T. At larger T, the decrease 
in Ao(T) compensates for this rise and eventually brings Iinei{T) to zero at Tc. Third, the 
behavior near Tc results from the fact that £'"'\T) oc Al{T) as Ao(T) goes to zero. 

The different forms for the critical currents discussed in this Section give rise to the 
possibility that measurements of jc(T) or the associated penetration depth may shed some 
light on the order parameter symmetry and the relative importance of the inter-layer hopping 
mechanisms in the cuprates. We have demonstrated that the effects of pairing symmetry in 
this simple model are quite strong at low temperatures: jc{T) oc T at low T if the pairing 
is (i-wave, and jc{T) is exponentially close to jc(0) if it is isotropic s-wave. This conclusion 
follows from the fact that j^''''^'^*(T) oc T and j™P(T) = j™*'^(T) = for (i-wave pairing, 
while all contributions are exponentially close to their T = values for isotropic s-wave 
pairing. Even in a more sophisticated model that includes anisotropy in the impurity- and 
boson- assisted scattering processes, the power-law dependence of jc{T) at low T for (i-wave 
pairing should persist due to the presence of the nodes in the gap function. Consequently, 
the low-temperature behavior of the c-axis critical current may be an important probe of the 
order parameter symmetry. In addition, we have seen that each inter-layer hopping process 
yields a jc with a distinct temperature dependence. Looking at the temperature dependence 
of the total jc{T) may therefore provide some clues about the relative importance of the 
different inter-layer hopping processes. The goal of the next section is to make a quantitative 
estimate of the critical currents for optimally doped and de-oxygenated YBCO and see if 
these expectations are realistic. 
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IV. RELATION TO EXPERIMENT 



The model we have presented builds on a natural connection between the normal- and 
superconducting-state properties. In this Section, we further exploit this connection by using 
experimental results on the normal-state properties of optimally doped and de-oxygenated 
YBCO to estimate jc{T) for these compounds. In the next two subsections, we describe the 
method for estimating the parameters in our model and extract the zero-temperature critical 
currents and penetration depths. We will use straightforward back-of-the-envelope reason- 
ing in making these estimates, but the magnitudes of these estimates are confirmed by more 
detailed fits to normal-state transport in this model.ilB In the succeeding subsection, we 
discuss the temperature dependence of the resulting critical currents and determine whether 
questions relating to the order parameter symmetry and the inter-layer hopping mecha- 
nism can be addressed in measurements of jc{T). Finally, we interpret recent experimental 
measurements based on the results of our calculations. 



Experimentally, the resistivity in optimally doped (Tc = 90 K) YBCO is roughly linear 
in temperature in both the in-plane and c-axis directions. Within the incoherent hopping 
model, we can attribute this behavior to purely direct hopping between the layers: from 
Eq. (^, the intra-layer scattering rate is directly reflected in the c-axis resistivity, implying 
that pab oc pc- In the superconducting state, this situation corresponds to jc(T) = j^^^'^^^T). 
To estimate the magnitude of this critical current, we observe that the band structure c-axis 
plasma frequency 



A. Optimally Doped YBCO 




(10) 



The T = limit of Eq. (|^) then yields 




(11) 
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Optical measurements give hujc = 400 meV,il from which we deduce jf""'\0) = 9 x lO'^ A 

_2 

cm . 

In layered materials, this critical current would also be related to the c-axis penetration 
depth measured in microwave experiments: Ac = yhcf/Sirecj^^^ , where c/ is the speed of 
light.0 This relation gives Ac = 0.5 yum for optimally doped YBCO, which is in reasonable 
accord with independent measurements of this quantity.ii In addition, Graf et al. estimate 
Ac = 0.9 /im for optimally doped YBCO, purely disorder-assisted hopping, and s-wave 
pairing,^ which is also in good agreement with our (independent) result. Assuming a density 
of states of 3 eV~^,0 we estimate t±_ ~ 10 meV. These results are summarized in the first 
column of Table |. 



B. De-Oxygenated YBCO 

To estimate the model parameters in the de-oxygenated case is more complicated, since 
we expect all inter-layer hopping processes to contribute to the critical current approximately 
equally. In particular, the boson-assisted hopping contribution should be significant in order 
to account for the upturn in pc at low temperatures. We begin by examining experimental 
measurements of the c-axis conductivity in de-oxygenated samples with = 40 K0 and 
find that 0"c ~ 1 ^l^^ cm^^ near Tc and dc is approximately linear in T at high temperatures 
with a slope ((icrc/(iT)highT ~ 0.018 cm~^ K^^. These asymptotic results are expected 
within our model, as can be seen from Eqs. (and see below). 

It is reasonable to assume (and detailed fits demonstrateSil) that the boson-assisted 
component of the normal-state conductivity has frozen out for T ~ T^, so we partition the 
remaining low-temperature conductivity equally between the direct and disorder-assisted 
processes. In addition, the lattice parameters for YBCoS give ctq = 0-24 fiil^^ cm~^ in 
Eqs. (0)-(^. Combining these facts with Eq. (^ immediately gives A^o^ ~ 1-4 x 10^'^ and 
j™P(0) = 3.7 X 10^ A cm~^. To extract the magnitude of the direct hopping, we note 
that optical measurements give h/rab = SvrAafe/csT with \ab = 0.2 to 0.4.0 Taking \ab = 
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0.3, one obtains fi/Tab = 6 meV at Tc = 40 K. Inserting this result into Eq. (0) produces 
A^o^i ~ 4 X 10~^ meV or 0.1 meV if iVo = 3 eV"^lll From NQt]_ and ctq, we obtain 

^■direct ^ 1.4 X 10^ A cm-2. 

At high temperatures, Eqs. (|^)-(^) can be written 

from which we obtain 

1 da, 



^ ' 4ec dT 



A2(0)/i„ei(T). (13) 

highT 



To compute this term, we also need the optical phonon frequency VIq. Recent c-axis-polarized 
Raman experiments indicate that the 500 cm~^ (720 K) 0(4) phonon is important in inter- 
layer hopping,!^ and so we take hQ^/ks = 720 K. Since H-Qq » Ao(0), Fig. |^ indicates that 
jinei^rp^ will monotonically decrease with T. Numerical computation of /inei(O) then gives 
f^^\0) = 12.5 X 10^ A cm-2. 

The estimated parameters for de-oxygenated YBCO are summarized in the second col- 
umn of Table |. Observe that t± has been reduced by two orders of magnitude from its fully 
oxygenated value: t± ~ 10 meV in = 90 K compounds, but t± ^ 0.1 meV in = 40 
K compounds. This strong variation in t± with doping is consistent with an exponentially 
decreasing wave-function overlap induced by the increasing c-axis unit cell dimension^ and 
by disorder and vacancies in the chains. Empirically, the c-axis resistivity increases by two 
orders of magnitude from optimally doped to de-oxygenated samples, 01311 which would also 
naively correspond to at least an order-of-magnitude decrease in t±. The data in Table | 
also indicate that jc"'^'(0) : f^^{0) : j^"'^^*{0) = 9:2.6:1, indicating that the Josephson current 
is dominated by the assisted processes in the de-oxygenated samples. This result is quite 
reasonable, given the importance of this term in determining the semiconducting tempera- 
ture dependence of Pc{T). As a final check to the consistency of our theory, we note that 
t± ~ Tc in optimally doped compounds, indicating weak incoherence, but t± « in the 
de-oxygenated compounds, which puts these compounds firmly in the incoherent regime. 
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C. Estimated jc{T) 



Using the parameters derived in the preceding two subsections in Eqs. (§)-(§), we can 
compute the total jc(T) in both optimally doped and de-oxygenated YBCO for both s- 
and d-wave paring. The results of this calculation are shown in Fig. ^. For the optimally 
doped case [Fig. |(a)], we see that the low-temperature jc(T) for d-wave pairing can be 
clearly distinguished from jc{T) for s-wave pairing, although jc(0) is independent of the 
pairing symmetry. In addition, the critical currents for either pairing symmetry are distinct 
from the conventional Ambegaokar-Baratoff result, and this difference should be observable. 
Specifically, the direct hopping dominates jc(T) in the fully oxygenated case, so the observed 
critical current should have the same temperature dependence as the in-plane superfluid 
density. Experimentally, this means that A^^(0)/A^j(T) should equal A^(0)/A^(T) in fully 
oxygenated YBCO. 

For the de-oxygenated case [Fig. §(b)], we again find a clear difference between the 
critical currents for s- and rf-wave pairing at low temperature (although it is not obvious 
from the figure). However, jc(0) now depends strongly on the pairing symmetry since the 
assisted processes are significant and contribute to the s-wave jc but not to the rf-wave 
jc- Additionally, while d-wave pairing yields a jc{T) very different from the Ambegaokar- 
Baratoff result, the s-wave critical current deviates from this form only very slightly. Thus, 
if the pairing is s-wave, it is unlikely that a measurement of jc{T) will enable one to deduce 
the relative importance of the inter-layer hopping mechanisms. 

The analysis of these figures allows us to draw two general conclusions. First, the order 
parameter symmetry should be indicated by the low-temperature behavior of the measured 
critical current regardless of oxygen content, being a power-law for (i-wave and exponen- 
tially close to a constant for isotropic s-wave. This result is robust and independent of the 
parameters estimated in the preceding two subsections. On the other hand, the reduction 
in magnitude of the critical current from s- to d-wave pairing in de-oxygenated samples 
depends sensitively on the model parameters, so it is unlikely that the pairing symmetry 
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will be indicated by an experimental measurement of jc(0), unless the observed magnitude 
is exceedingly small (indicating d-wave pairing). Second, the presence of incoherent hopping 
in Ac should be revealed by the observation of either a 3D BCS character if the pairing 
is d-wave or a 3D BCS character that crosses over as a function of hole doping to a 2D 
Ambegaokar-Baratoff character if the pairing is s-wave. It is unlikely, however, that these 
measurements will in isolation allow one to disentangle the relative contributions of the 
different inter-layer hopping processes in the de-oxygenated samples. 



D. Comparison to Experiment 

Our model can be compared to direct experimental measurements of jc{T) or to indirect 
measurements of this quantity through, for example, microwave penetration depth measure- 
ments. In the latter experiments, the penetration depth Ac is related to the critical current 
jc by the relation jc{T)/jc{0) = A^(0) / A^(T).i§ Given the results of the preceding subsection, 
we should ideally compare our theoretical estimates to measurements of the critical current 
or penetration depth in YBCO. Unfortunately, the experimental situation is not yet resolved 
in this cuprate: some penetration depth measurements show an SNS-like A^(0)/A^(T) which 
is incompatible with our model,0 while others yield more conventional curves.^ However, 
we may still examine measurements on other cuprates and look for some of the generic 
features of our model. 

The most thoroughly studied cuprate with regard to intrinsic Josephson junction ef- 
fects is BiaSraCaCuaOg (BSCCO).0 Direct critical current measurements are available for 
this material that show Ambegaokar-Baratoff behavior in oxygen-annealed samples and a 
temperature-independent jc{T) at low T that falls rapidly to zero near Tc in argon-annealed 
samples.0 Within our model, the Ambegaokar-Baratoff-like jdT) in the oxygen-annealed 
samples suggests s-wave pairing and a significant contribution to jc from assisted hopping, 
which is analogous to de-oxygenated YBCO [cf. Fig. ^(b)]. The approximately constant 
jc{T) in the argon- annealed samples is also indicative of s-wave pairing, but the rapid re- 
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duction of jc near Tc is difficult to explain within our model. 

In addition to BSCCO, recent measurements of the c-axis penetration depth on 
La2-xSrxCu04 (LSCO) have been reported and interpreted within an intrinsic Joseph- 
son effect picture. As with oxygen-annealed BSCCO, LSCO shows an approximately 
Ambegaokar-Baratoff jc{T), which we interpret as evidence for significant assisted inter- 
layer hopping as well as s-wave pairing. However, the data do not extend to very low 
temperatures, so additional experimental measurements in this compound (and in BSCCO) 
would be illuminating. 

V. CONCLUSIONS 

In this paper, we have considered an incoherent hopping model that reproduces the 
features of the c-axis resistivity in both optimally doped and underdoped cuprates. We 
have assumed that there are three processes which transport quasiparticles from one Cu02 
layer to another: a direct hopping induced by wave-function overlap, an impurity-assisted 
hopping due to disorder, and a boson-assisted hopping induced by, for example, modulation 
of the inter-layer distance by phonons. By performing calculations to second order in the 
inter-layer hopping amplitudes, we have shown that this model yields an intrinsic Josephson 
effect in the superconducting state, and we have analyzed the contributions to the Josephson 
critical current arising from these three processes. We find a robust dependence of the 
low-temperature c-axis critical current on the order parameter symmetry which should be 
reflected in experimental measurements oi jc{T) or the c-axis microwave penetration depth 
ratio Xm/^l{T). 

This conclusion is supported by the application of our model to YBCO, where we have 
estimated the magnitude and temperature dependence of jc and find a clear-cut signature 
of the pairing symmetry in the low-temperature behavior of jc{T) for both fully oxygenated 
and de-oxygenated compounds. Additionally, our model predicts that the c-axis penetration 
depth ratio should have either a 3D BCS character if the pairing is d-wave or a 3D BCS 
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character which crosses over as a function of hole doping to a 2D Ambegaokar-Baratoff 
character if the pairing is s-wave. Measurements on BSCCO and LSCO suggest that the 
pairing in these materials is s-wave and the inter-layer hopping has a significant assisted 
component, but direct comparison of our estimates with YBCO is not possible because the 
data have not yet converged. Further experimental work on both fully oxygenated and 
de-oxygenated YBCO is clearly required and would provide a stringent test of our theory. 
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FIGURES 

FIG. 1. Feynman diagrams used to compute the c-axis conductivity and Josephson critical 
current within the incoherent hopping model. The solid dots at the vertices represent the direct in- 
ter-layer hopping amplitude t± , the dashed line with a cross represents the averaged inter-layer dis- 
order, and the wavy line represents the propagator for boson-assisted hopping. The two electronic 
Green's functions in each diagram correspond to nearest-neighbor layers; normal-state propagators 
are used to compute the conductivity while Gor'kov propagators are used in the superconducting 
state. See text for details. 

FIG. 2. Normalized contributions to the c-axis critical current jc{T)/jc{0) as a function of 
the reduced temperature T/T^ from various inter-layer hopping processes and pairing symmetries. 
Shown are critical currents for d-wave pairing and direct hopping (j^"'*''^*(T)/j^"'°'^*(0) from Eq. (|^, 
dashed line) , s-wave pairing and direct hopping (jf '^'*(T)/jf '■^^*(0) from Eq. (|), dot-dashed hue), 
and s-wave pairing and impurity-assisted hopping (jc™^(^)/ic™^(0) from Eq. (^), solid line). These 
results are identical to, respectively, the in-plane BCS superfluid density for d-wave (dashed line) 
and s-wave pairing (dot-dashed line) and the critical current in a macroscopic SIS tunnel junction 
with s-wave pairing obtained within the Ambegaokar-Baratoff (AB) theorylll (solid line). Note 
that the d-wave impurity-assisted hopping critical current vanishes for isotropic scattering. 

FIG. 3. Normalized boson-assisted c-axis critical current /inei(^) [Fq. (^] as a function of 
temperature T for Einstein phonons with frequencies 7iJ7o/^B (from top to bottom) of 800 K, 
400 K, 200 K, 100 K, 50 K, and 25 K with fixed coupling constant A (solid lines; see text for 
details). The temperature dependence of the gap function used in this calculation is determined 
from the s-wave BCS equation with a Debye temperature of 1000 K and a coupling constant fixed 
so that the critical temperature is 100 K. These choices yield Ao(0) = 178 K. For comparison, the 
Ambegaokar-Baratoff (AB) result for an SIS Josephson junctionEl normalized to the largest /inei(O) 
is also shown (dot-dashed line). 
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FIG. 4. Estimated c-axis Josephson critical current jc{T) obtained within the incoherent 
hopping model as a function of temperature T for (a) optimally doped (Tc = 90 K) and (b) 
de-oxygenated (Tc = 40 K) YBa2Cu307_5 with s-wave (dot-dashed line) and d-wave (dashed line) 
pairing. For comparison, the Ambegaokar-Baratoff resunjll normalized to the largest jc(0) is shown 
as a solid line. The parameters used to generate these curves are discussed in the text. 
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TABLES 

TABLE I. Estimates of the direct hopping amphtude t±, the zero-temperature critical currents, 

and the zcro-tcmpcraturc penetration depths for fuUy oxygenated (Tc = 90 K) and de-oxygenated 
{Tc = 40 K) YBa2Cu307_5 within the incoherent hopping model. The critical currents are pre- 
sented for the direct {jc"^^^{0)), impurity-assisted (jc™^(0)), and boson-assisted {j™^^{0)) inter- layer 
hopping processes along with the total critical currents (j^(0), Jc (0)) ^^'^ corresponding penetration 
depths (A^(0), (0)) for s- and d-wave pairing. See text for details. 



Optimally Doped 



De-oxygenated 



t± (meV) 



^•direct ^0) (A Cm-2) 



jrP(O) (A cm-2) 



10 

9 X 10^ 
~ 

~ 



0.1 
1.4 X 10^ 
3.7 X 10^ 
12.5 X 10"^ 



jm (A cm-2) 
j^(0) (A cm-2) 



9 X 10^ 
9 X 10^ 



17.6 X 10^ 

1.4 X 10^ 



A?(0) (/xm) 



0.5 
0.5 



11 
40 
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